function dy = ressler(t, y)
dy = zeros(3, 1);
dy(1) = -y(2) - y(3);
dy(2) = y(1) + 0.2 * y(2);
dy(3) = 0.2 + y(3) * (y(1) - 5.7);
end

